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Abstract: Differential rates in the decay B® — > J/ij)(f> with (f> — > K + K~ and J/if) — > fi + fi~ 
are sensitive to the CP-violation phase j3 s = arg ((—Vt s V t * b ) / (V cs V* b )) , predicted to be very 
small in the standard model. The analysis of — > J/ijjcj) decays is also suitable for 
measuring the B® lifetime, the decay width difference Ar s between the B® mass eigenstates, 
and the B® oscillation frequency Am even if appreciable CP violation does not occur. In 
this paper we present normalized probability densities useful in maximum likelihood fits, 
extended to allow for S- wave K + K~ contributions on one hand and for direct CP violation 
on the other. Our treatment of the S"-wave contributions includes the strong variation of 
the S-wave/P-wave amplitude ratio with m{K + K~) across the (j> resonance, which was 
not considered in previous work. We include a scheme for re-normalizing the probability 
densities after detector sculpting of the angular distributions of the final state particles, and 
conclude with an examination of the symmetries of the rate formulae, with and without 
an S-wave K + K~ contribution. All results are obtained with the use of a new compact 
formalism describing the differential decay rate of B® mesons into J/ipcfi final states. 
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1. Introduction 

The decay of B® — > J/ip<p, a transition of a pseudoscalar into two vector mesons can be 
thought of as six independent decays. The initial B® system consists of a heavy and a light 
mass eigenstate, and the J/ip(j) system to which it decays is characterized by three distinct 
orbital angular momentum states. A maximum amount of information about this system 
can be obtained from analyses which disentangle the two initial states and the three final 
states. The experimental technique of flavor tagging infers a meson's flavor at production 
time as B® or B® and is the key to disentangling the two initial states. Flavor-tagged 
B® — > J/tp(p decays are of great interest in particle physics because of their sensitivity to 
the CKM phases [1] and to anomalous mixing phases from physics beyond the standard 
model [2]. Recently the CDF and DO collaborations have constrained the CKM phases 
in both untagged analyses [3, 4], and flavor-tagged analyses [5, 6]. These analyses are 
based on complete differential rates for the decay given in Ref. [7]. They use the angular 
distributions of the decay products to disentangle the three final states. 

In this paper we re-express the differential decay rates in Ref. [7] , using a new formalism 
that makes explicit a number of symmetries that are otherwise hidden. These formulae are 
then extended to the case in which the final state in the decay B® — > J/ip(f> includes decays 



of type Bg — > J/ipK + K~ (kaons in an S- wave state), which has been suggested [8] to be 
an important effect. After including the S'-wave contribution in the theoretical description, 
we identify the symmetries of the modified formulae. 

In addition to S'-wave effects, we also investigate other aspects of the differential decay 
rate formulae. We include the effects of possible direct CP violation. In addition we show 
how interference between CP odd and CP even J/ip<ft final states effectively tags the flavor 
of the Bg meson at decay, allowing for the possibility to observe B® — > B® flavor oscillations 
in a flavor-tagged analysis, even in the absence of CP violation effects. 

Experimentally, the differential rate formulae are used to construct likelihood func- 
tions based on normalized probability density functions (PDFs). In this paper we include 
normalization constants where appropriate in all expressions for transition amplitudes and 
PDFs. Detector angular acceptance is an important effect which must be included in these 
probability densities. However, the inclusion of this effect disturbs the normalization of the 
PDF. We present a scheme for normalizing the probability density analytically, as required 
for unbinned maximum likelihood fits. 

2. Phenomenology of the B Q S — > J/'tp(fi Decay 

We first summarize the phenomenology of the B® system and the decay B® — > J/tp'P ~~ 
/i + /i~ K + K~ . Two flavor eigenstates, \B®) and \B^), mix via the weak interaction. The 
two mass eigenstates 

\Bf) = p\B° s ) - q \B^) = p \B° S ) + q \B° S ) 

are labeled "heavy" and "light". The mass and lifetime differences between the and 
Bg states can be defined as 

Am = m H - m L , AT = T L - F H , T = (Th + T L )/2 , 

where mn,L and Th,l denote the mass and decay width of Bf and B^ (with this definition 
both Am and Ar are expected to be positive quantities). The heavy state decays with a 
longer lifetime, th = 1/Ti7i while the light state decays with the shorter lifetime tl = l/r^, 
in analogy to the neutral kaon system. The mean lifetime is defined to be r = 1/T. 
Theoretical estimates predict AT/T to be on the order of ~ 15% [2]. Linear polarization 
eigenstates of the J/ijj and (j) provide a convenient basis for the analysis of the decay [9]. 
The two vector mesons can have their spins transversely polarized with respect to their 
momentum and be either parallel or perpendicular to each other. Alternatively, they can 
both be longitudinally polarized. We denote these states as \P\\), \P±) } and \Vq). 

In the standard model, CP violation occurs through complex phases in the CKM 
matrix [10]. Large phases occur in the matrix elements V u b and Vtd- While these matrix 
elements generate large CP violation in the B° system, they do not appear in leading order 
diagrams contributing to either B® f-> B® mixing or to the decay B® — > J/ipcp. For this 
reason the standard model expectation of CP violation in B® — > J/ipcp is small. In the 
limit of vanishing CP violation, the heavy, long-lived mass eigenstate B^ is CP odd and 
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decays to the CP-odd, L=l orbital angular momentum state \V±). The light, short-lived 
mass eigenstate Bjf is CP even and decays to both CP-even L=0 and L=2 orbital angular 
momentum states, which are linear combinations of |Po) an d |P||)- 

The small CP violation in B® — > J/ip4> can be quantified in the following way: we 
define Ai as the decay amplitude {B s \H\Pi) and Ai as the decay amplitude {B s \H\Pi) 
where i is one of {||,_L,0}. All CP observables in the system are characterized by three 
quantities Aj = In the standard model the Aj are given as Ai = ±exp(i2/3 s ) where 

the positive and negative sign applies to the CP even and odd final state, and 

The standard model expectation [II] is 2f3 s = 0.037 ± 0.002, a very small phase which does 
not lead to appreciable levels of CP violation. New physics can alter the mixing phase, 
while leaving A very nearly unimodular. In this paper we consider, however, also the case 
in which |A| 7^ 1. 

3. Differential Rates 

The state of an initially pure B® or B® meson after a proper time t has elapsed is denoted 
as \B® phys (t)) and \B^ phys (t)) . Transitions of these states to the detectable n + fi~K + K~ 
can be written as 

(MV K + K~\H\Bl phys {t)) 

= {^^K + K~\H\Vi) (Vi\H\B° s ) (B° s \Bl phys (t)) 

i 

+ Y,(V + V~ K+K ~\ H \K) {Vi\H\m) (B° s \B s>phys (t)), 

i 

(MV K +K-\H\B° Sjphys (t)) 

= Y.i^^K+K-lHlV.) (Vi\H\B° s ) (B° s \B° s>phys (t)) 

i 

+ Y,^ + ^ K+K ~\ H \^) <Pi\H\m) (B° s \Bs, P hys(t))- 

(3.1) 

where H is the weak interaction Hamiltonian. The expression can be written much more 
simply, by defining time-dependent amplitudes for |P[?) and \Bg) to reach the states \Pi) 
either with or without mixing: 

Ait) = (Vi\H\B° s ) (B° s \Bl phys (t)) + {Vi\H\B2) (B° s \Bl phys (t)) , 
Mt) = (Vi\H\B s ) (B° s \B° S)phys (t)) + (Vmm) (B° s \B SjPhys (t)) ■ 

Then: 

(^^K + K~\H\Bl phys (t)) = ^2Mt)e~ imt {^^K+K^HWi) , 
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{^^-K+K-\H\Bl phys {t)) = Y<MtY- imt (n + H-K + K-\H\Vi 



(3.2) 



where the time dependence of Ai(t) and Ai(t) is: 

e -rt/2 

Ai(t) 



\/t h + t l ± cos 2(3 s (t l - t h ) 

e~ rt / 2 

y/r H + T L ± cos 2f3 s {t~ l - t h ) 



E + {t)±e 2i ^E^{t) 
±E + (t) + e~ 2ilSa E-{t) 



(3.3) 



and where the upper sign indicates a CP even final state, the lower sign indicates a CP odd 
final state, 



E±(t) = 



l r 



,+( 



0* 



' -Ar . ■ Am % 



and the Oj are complex amplitude parameters satisfying: 



Ski 



i. 



(3.4) 



(3.5) 



The final state K + K~ is characterized by three decay angles, described in a 

coordinate system 1 called the transversity basis [1]. In the J ftp rest frame, the x-axis is 
taken to lie along the momentum of the <p and the z-axis perpendicular to the decay plane 
of the <p. The variables (6, <p) are the polar and azimuthal angles of the fi + momentum in 
this basis. We also define the angle tp to be the "helicity" angle in the <p decay, i.e. the 
angle between the K + direction and the x-axis in the <p rest frame. With these definitions, 
the muon momentum direction in the J/tp rest frame is given by the unit vector 



h = (sin 9 cos <p, sin 8 sin ip, cos 9) . 
Let A(i) and A(t) be complex vector functions of time defined as 

A(t) = Ao(t) cos ip, ll ^A= ,t 



(3.6) 



V2 



V2 



Mt) 



t u . , ^||WsinV .^ ± (t)sin^ 
Ao(t) cosip, t= , %- 



V2 



V2 



(3.7) 



where Ai{t) have now been normalized. For experimental measurements we are concerned 
with normalized probability density functions Pb and Pg for B and B mesons in the 
variables t, cosip, cos 6, and tp, which can be obtained by squaring Eq. (3.2). The formulae 
of Ref. [7] are then equivalent to: 

9 



P B (P,<p,il>,t) 
P s (8,tp,*p,t) 



-\A(t) x n| 2 
16tt' w 1 

— |A(t) x n| 2 
16vr' W 1 



(3.8) 



1 An alternate basis called the helicity basis is discussed further in Section 9. 
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which give a picture of a time-dependent polarization analyzed in the decay 2 . The factors 
of 9/167T are normalization constants, and are present in order that 

/ Pj(ip,d,V,t)d(cosi>)d(cos0)dipdt = 1. (3.9) 

j=B,B 

The quantities |aj| 2 give the time-integrated rate to each of the polarization states. The 
values of Ai(t) at t = will be denoted as Ai. To translate between the o's and the A 7 s 
one can use the following two sets of transformations: 

14 |2 \ n |2 \ A ±\ 2 

\ A ±\ = i , / m — 19 I°-lI "~ 



l + (y-l)\a ± \ 2 

= i + ( y -i)K| 2 |a »' 

i, ,2 _ l Q o! 2 I |2 

A ° ~~ TT7 TYi 12 a ° 

1 + (y- 1) a± 2 



y + (i- 


y)|A ± | 2 




l 2 y 


y + (i- 


y)l^±l 2 




l 2 y 


y + (i- 


y)l^4±l 2 



(3.10) 



where y = (1 — z)/(l + z) and z = cos 2/3 s Ar/(2r). The relation (3.5) insures that 

i 

Eq. (3.8), together with the definitions in Eqs. (3.3), (3.4), and (3.6) can be used as a decay 
model for an event generator, and is suitable for use as a fitting function in the absence of 
detector effects. 

4. Detector Efficiency and Normalization 

The detector efficiency e(tp,6,ip), when introduced into the above expression, disturbs the 
normalization of Eq. (3.9). We restore it by dividing by a normalization factor N, 

P '(rp, 9, ip, t) = ip(^, d, ip, i)ety>, 6, if) , 

N = j ^2 Pi(ip,d,if,t)e(tp,e,ip)d(cosip)d(cos9)dcpdt. 

J i=B,B 

(4.1) 

Suppose that the efficiency e(ip } 6, ip) can be parametrized as 

ey>,0,<p) =cf m P k (co S ^)Y lm {e,<p), (4.2) 

where cf m are expansion coefficients, P k (cos ip) are Legendre polynomials, and Yi m (0,<p) 
are real harmonics related to the spherical harmonics through the following relations: 

Yim = Y™ (m = 0) , 

Y lm = ±(Yr + (-l) m Yf m ) (m>0), 

Y ^ = ^l™ 1 ~ (-l) |m| iT H ) (m < 0) . (4.3) 



2 Throughout this paper, when writing the dot product of two complex vectors, we always imply complex 
conjugation on the second operand. 
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The products Pfe(cos ip)Yi m (9, <p) constitute an orthonormal basis for functions of the three 
angles. The detector efficiency (obtained, for example, from Monte Carlo simulation) can 
be fit to the first few of these polynomials. A straight-forward calculation shows that: 



N 



3 

8J^ 



4r° 
4C 00 



(kol 



Ar 2 



+ 



15 

3 



hil + W±\ ) 
-Kl 2 ) 



8V5vr 



^(|a | 2 + |an| 2 -2|a ± | 



+ 



4c 



20 



15 



+ KI 2 ) 



1, 

2 11 

sin2/3 s (TL - t h ) 



IQV^M y/(( TL - t h ) sin 2/3 s ) 2 + At l t h 



(a|a_L + a||ol)(-c^_ ] 



15 



c 2 -iJ 



+ 



9 V2 



sin2/3 s (r L - t#) 



16 ^15^ ^(( TL - t h ) sin 2/3 s ) 2 + At l t h 



(a^a± + a a* ± )( 



32 



+ ... 



+ 
+ 



9 



8^/15^ 
9 V2 



2d 



4 r 2 

22/ I 1 2 ) I |2\ 4C 22/|„ |2 , ^ |2> 

\—\ a o\ + a|| J 



1 



lT(N J + 2i a ii 



7T 



(aoay + a a|)( 



, 7TC. 



2-2 



ttc: 



2-2 



8 



32 



+ ... 



(4.4) 



The numerical factors +7r/8 and — vr/32, appearing together with and C2 _2 in the 
infinite series, are the integrals 



Pfc(cosV0 cos(V') sin ipd(cos ip) 



(4.5) 



While this series is infinite, the number of basis functions needed to fit detector efficiencies 
in a particular analysis is finite and determined chiefly by the size of the data sample. With 
the factors in Eq. (4.5) the normalizing factor can be adapted to account for all terms used 
in the expansion of the efficiency. Eq. (4.4) represents an analytic normalization of the 
fitting function and provides an efficient way to compute the likelihood during a maximum 
log likelihood fit. The orthonormality of the basis functions has been used to reduce the 
expression to its final form. 



5. Time Development 

The short oscillation length of the B® meson [12], 2irc/Am ~ 106 fim, requires us to 
account for resolution effects when fitting the rates of flavor-tagged decays, even using the 
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best silicon vertex detectors, which have proper decay length resolutions on the order of 
25 fim. Certain time-dependent functions arising from particle-antiparticle oscillations, 
particularly those expressed as the product of exponential decays and harmonic functions 
with frequency Am, must be convolved with one or more Gaussian components describing 
detector resolution. This convolution can be carried out analytically, using the method 
described in Ref. [13] for the evaluation of certain integrals which are equivalent to complex 
error functions. In this step one requires that various components of the time dependence 
first be separated from Eq. (3.8). The time development of Ao(t) and A\\(t) amplitudes 
are identical, but differs from that of A±(t). We begin by decomposing 



where 



A(t) = A+(t) + A_(t), A(t) = A+(f) + A_(t) 

ay sin ijj 



A+(i) = A+/+(t) = (aocosV>, 



V2 



.0) •/+(*) 



and 



_ an sin'0 

A+ (t) = A+/+ (t) = (a cos if), - -L=- , 0) • /+ (t) 



A_(t) = A_/_(t) = (0,0, i^^-) ■ /_(t) , 
A_(t) = A_/_(t) = (0,0, • , 



and we define 

f±(t) 

f±(t) 



-rt/2 



+ T i ± COS Ws (jL ~ T H ) 

e- Tt l 2 

y/r H + T L ± cos 2f3 s (t l - t h ) 



V2 

E+{t)±e 2i ^E^{t) 



We then have in place of Eq. (3.8) 

{|A+(t) x h\ 2 + |A_(t) x n| 2 + 2ifes((A+(t) x n) • (A* (t) x n))} 



16vr 
9 

16tt 



{|A + xn| 2 |/+W| 2 + |A-xn| 2 |/-(t)| 2 
+ 2Re((A + x n) • (A* x n) ■ f + (t) ■ /*(t))} 



and 



P 5 (6,i/>,<p,t) 
9 



16vr 
9 

167T 



(|A+(t) x n| 2 + |A_(t) x n\ 2 + 2i?e(A+(f) x n) ■ (A* (t) x n))} 
{|A + xn| 2 |/ + (t)| 2 + |A_xn| 2 |/_(t)| 2 



+ 2i?e((A + x n) • (A*_ x n) • f + (t) ■ f*_(t)} 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



(5.6) 
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where (for B) the diagonal term in Eq. (5.6) is 

- 2 _ l(l±cos2/3 s )e- rL * + (l=Fcos2/3 s )e- r «*±2sin2/3 s e- r *sinAmi 
1 Wl ~2 r L (l ± cos 2(3 S ) + t h (1 t cos 2f3 s ) ' 

(5.7) 

while (for B) the diagonal term in Eq. (5.5) is 

2 _ 1 (1± cos 2/3 s )e" r ^ + (1 =F cos 2/3 s )e- r »*=F 2 sin 2^ s e~ rt sin Amt 
1 Wl ~2 r L ( 1 ± cos 2/3 s ) + t h ( 1 T cos 2/3 s ) 

(5.8) 

and (for £?) the cross-term, or interference term in Eq. (5.6) is 

-,.-», -e" rt cos Amt - icos2f3 s e~ vt sin Ami + zsin2/3 s (e" ri * - e~ rH *)/2 

/+(*)/- W — / > 

y [(tl - th) sin 2/3 s ] 2 + 4t l t^ 

(5.9) 



while (for B) the interference term in Eq. (5.5) is 

e~ rt cos Amt + i cos 2/3 s e~ rt sin Amt + z sin 2(3 s (e~ FLt - e~ Vllt )/2 



U(t)f-(t) 



[(tl ~ t h ) sin 2/3 s ] 2 + At l t h 



(5.10) 



This accomplishes the desired separation. In the fitting function, to accommodate the 
proper time resolution, one has only to replace all time-dependent functions with their 
smeared equivalents. 

6. Sensitivity to Am 

It can be noticed that the time development of the interference term, expressions 5.9 and 
5.10, contain undiluted mixing asymmetries even in the case of no CP violation, i.e., when 
P s = 0. Let us try to better understand the mechanism by which the flavor of the B® 
meson is tagged at decay time, by first rewriting Eq. (3.1) using the Bg and Bg states in 
the expansion rather than the B® and B® states: 

( M +/i- K+K-\H\Bl phys (t)) = 

^ (n + ti-K + K-\H\Vi) {Vi\H\Bf) (Bf \B® phys {t)) 

i 

+ ]T (ji+vL-K+K-lHlPi) (Vi\H\B^) (B^\B° s>phys (t))- (6-1) 

i 

Now, we take the limit of zero CP violation in the B® system, such that {V\\\H\Bg ) = 
(Vo\H\Bg ) = (V±\H\Bg) = 0, and only three of the six terms in Eq. (6.1) remain: 

K+K-\H\Bl phys (t)) = 

(fi + ^K + K-\H\V ± ) (V±\H\B?) (B?\Bl phys (t)) 
+ (n + ^K + K~\H\V ) (Po\H\Bt) (B^\B° S)phys (t)) 
+ {ja+v-K+K-IHIPu) {P\\\H\Bf) (B^\B° Stphys (t))- (6-2) 
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When the expression is squared, the interference terms are the cross terms involving both 
the product of a CP-even and a CP-odd amplitudes. The time dependence of these terms 
is contained in the factor: 



j [{Bf\B° SjPhys {t)) + {Bg\B G sphys {t))) 2 
- {{ B f \B° Sjphys {t)) - (Bs\B® phys (t))) 




V2 



B 



,0 



) 



21 



(Bs\Bg phys (t)) — {Bg\Bg phys (t)) 



(6.3) 



This factor takes the value +1/2 when the meson is pure B®, and -1/2 when the meson 
is pure B®, and in general oscillates between these two values. Thus the interference term 
effectively tags the flavor of the B® at decay. This provides a way to observe B® — > B® 
flavor oscillations using a sample of flavor-tagged B® — > J/ip(f) decays which can be collected 
with a simple dimuon trigger. This may open a particularly interesting avenue for the LHC 
experiments to observe B® mixing using a J/tp trigger. 

7. Incorporating Direct CP Violation 

An asymmetry either in the decay rate (|j4j/Aj| ^ 1) or in the mixing (|^/p| 7^ 1) such 
that |A| 7^ 1 is direct CP violation. In the case of direct CP violation A does not lie on 
the unit circle in the complex plane, and we need two parameters to describe it which we 
will take to be C = Re(X) and S = Im(X). Experimentally, even if one sets out to extract 
(3 S assuming the constraint |A| = 1, it is nonetheless of interest to test that constraint, 
since sensitivity to C and S arise from very different features of the detector. In that case 
we must revisit not only the functional form of the differential decay rates, but also the 
normalization. The amplitudes in Eq. (3.3) must now be written as: 



Ai = M±e- rt/2 [E+(t) ± XE- (t)] at 
Ai = N±e- Vt / 2 [±E+(t) + E-(t)/\] 



'i j 



(7.1) 



where 



A4 




+ 



1 + Am 2 r 2 



T 



±45- (1-|A| 2 ) Arm- -2(1- |A| 2 ) 2 



} 



2 



- 9 - 



These amplitudes can readily be seen to reduce to those of Eq. (3.3) in the limit of C 2 ±5 2 = 
|A| 2 — > 1. The normalization of detector efficiency, Eq. (4.4), becomes: 



N 



8^ 



%|a | 2 + h| 2 + K| 2 ) 



+ 4 | , (2|ao| 2 



+ 



8V57T 

Ac 2 



2c! 



(|a | 2 + |a||| 2 -2|a ± | 2 ) 



1 



5 (\ao\ 2 -- 2 H 2 + \ ai _\ 2 ) 



9 



f6V15vr 



N+NS ■ (t l - TH ) 



(a|a_L + a||al)(-c°_ 



15 



c 2 -iJ 



V -^^Ln + NS-(t l -t h ) 



16 715^ 



(aga_L + a a* ± )( 



TIC 



21 



7rc: 



21 



32 



+ ...) 



+ 



9 



zc 22/ I 1 2 



8\/l57r 
9 V2 
16 v/l5^r 



(-|o r + hil) 



4c 



22 



15 



(|ao| 2 + da,,| 2 ) 



Ki a || +a a|)( 



1 3 



32 



+ ... 



Finally, the explicit time development, Eqs. (5.7), (5.8), (5.9) and (5.10), must be replaced 
with the more general forms: 



|/±(<)| a 
\f±(t)\ 2 

/+(*)/-•(*) 
/+(*)£(<) 



ATI 
4|A| 2 



[((1 + |A| 2 ) ± 2C)e- VLt + ((1 + |A| 2 ) T 2C)e^ Ht 



+ (±45 sin Ami - 2(1 - |A| 2 )cosAmi) e 



-rti 



AT: 



± 



[-e~ vt (2(1 + | A| 2 ) cos Ami + 4iCsinAmt) 



[((1 + |A| 2 ) ± 2C)e- TLt + ((1 + |A| 2 ) =f 2C) e " r ^ 

- (±45 sin Ami - 2(1 - |A| 2 )cosAmt) e~ Tt ] , 
M+M - 

+e -?Lt ( (1 _ | A |2) + 2iS j + e -T H t ( (1 _ | A | 2) _ 2 -^j j 

A/ "+^ Vl [e" r * (2(1 ± |A| 2 )cosAmt±4iCsinAmt) 
+e~ TLt ((1 - |A| 2 ) + 2iS) + e~ VHt ((1 - |A| 2 ) - 2iS)\ , 



which can be seen to reduce to expression 5.7, 5.8 and 5.9, 5.10 as |A| 2 — > 1. 
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8. Incorporating a Contribution from B® — y J/ipK + K (Kaons in an S- 
Wave State) 

It has been suggested [8] that a contribution from ,5- wave K + K~ under the (ft peak in 
B® —y J/ift(ft decay may contribute up to 5-10% of the total rate. A normalized probability 
density for the decay B® — y J /iftK + K~ (kaons in an S- wave state) can be worked out by 
considering the polarization vector of the J ftp in the decay and proceeding as in [9]. The 
resulting expressions 



Q B (0,(f,ip,t) 



16tt 
3 

16tt 



B(t) x n| 



B(t) x n 



- |2 



do not depend at all on the angle ift (which is the helicity angle in the (ft decay). In the 
previous expression 



B(t) = (B(t), 0,0) , 
B(t) = (B(t), 0,0) 



(8.2) 



where the time-dependent amplitudes, 

e -rt/2 



B(t) 
B(t) 



\Jt~h + t l - cos 2f3 s (t l - T H ) 

e- Tt l 2 

\/th + t l - cos 2f3 s (t l - r H ) 



E + (t) 



E-(t) 



-E + (t) + e~ 2if3s E^(t) 



i.3) 



reflect the CP-odd nature of the J/ipKK final state. 

When both P-wave and 5-wave are present, the amplitudes must be summed and then 
squared. The P wave has a resonant structure due to the ^-propagator, while the S-wave 
amplitude is flat (but can have any phase with respect the P-wave). Suppose that in our 
experiment we accept events for which the reconstructed mass m(K + K~) = \i lies within 
a window [L\ < fi < fifa. The normalized probability in this case is 



p B (9,(f,^,t,n) 



9 
16vr 
9 

16tT 



y/l-F a g{ t i)A.{t) + e iS >yfF l 



V3 



B(t) 



V3 



x n 



x n 



where we use a nonrelativistic Breit-Wigner to model the (ft resonance 3 

1 



TJ2 



Au /j, \x§ + iT (j> /2 



(8.4) 



(8.5) 



3 We shall have more to say about that, later. 
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a flat model for the S-wave mass distribution 

MM) = 7/|= (8-6) 

and define 

, _i 2(/i/„ - fj,<f>) , _i 2 (Wo - /q y \ 

u hi = tan taio = tan (8.7) 

and 

A/i = n h i - m Aw = id hi - ujio . (8.8) 

In these equations, F s is the S-wave fraction; /j,^ is the (p mass (1019 MeV/c 2 ); Fj, is the cf> 
width (4.26 MeV/c 2 ), and 5 S is the phase of the S'-wave component relative to the P-wave 
component. 

In the presence of an S'-wave contribution, the normalization of Eq. (4.4) must be 
generalized; in order to do this we first define the quantities 



F s (l - F s )r e~ i5s 
^ = V 2A^Aw ' + ziy 2 (8 ' 9) 



and 



' f 3 (i - F s )r _ ^_ iSs _ 1q /ifei - ^ + ^W2 



(8.10) 



Then the normalizing factor appropriate for Eq. (8.4) is 

M = (1 - F s ) ■ N + 2i?e [X M • iV'] + F s ■ N" (8.11) 
where iV is given in Eq. (4.4), and 

N' = v / 3*a^(— + — A=4o" — ^ c 2 2 ) 

DV 71 ^ 12V57T 4V157T 

3/2 at ,7T n 7T 9 

+ 16V^ ( 2 C - 2 -8 C - 2 + - ) 

16V 5^ ± ^((r L -T H ) S m2f3 s y + ir L T H y 2 21 8 21 



.12) 



N" = ~c° + -^=c° - 3 4, . (8.13) 



and 



The numerical factors +ir/2 and —ir/8 appearing together with C2 1 and _ 2 in the infinite 
series are the integrals 

Pk(cosi/j) sin?/>d(cos ip) . (8-14) 
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We now work out the explicit time and mass dependence of the differential rates. We 
will use Eq. (5.5) together with the analogous equation for the pure S- wave differential 
rate: 



Q B (0,il>,<p,t) = y^IBW x "I 2 



Jt|Bxn| 2 |/_(t)| 2 (8.15) 



and 



Q S {W,<P,t) = y|^|B(£) x n\ 2 

= J-|B x h\ 2 \Ut)f . (8.16) 

where the vector B = x = (1, 0, 0). The full probability densities, which can be used in a 
time-, angle-, and <p mass-dependent fit, are obtained by expanding Eq. (8.4). We get 

PB (0,Tp,if,t,fl) = 

(1 -F.&£- * ■ P B {B, V, <p, t) 

+F a —Q B (e,il;,<p,t) 



+2j^Re ((A_ x h) ■ (B x h) ■ \f-(t)\ 2 

+(A+ x n) • (B x n) •/+(*)•£(<))] 



.17) 



and 



Pb (d,ip,ip,t,n) = 

(1 - F.) ^ • - -i , • P B (9, V», t) 

+Fs-£-QBM,<P,t) 



\/27 

+2j 6n Re ((A_ x n) • (B x n) • |/_(i)| 2 

+(A + xn)-(Bx ft) •/+(*) ■/*(*))] . 



(8.18) 



- 13 - 



In case one does not want to observe the 0-mass variable fi, one can integrate it out. 
Then one obtains 



PB (0,ip,<p,t) = 

(1 - F.) ■ P B (6, A if, t) + F S Q B (6, <p, t) 

+2^i?e [X M ((A_ x h) • (B x n) ■ |/_(t)| 2 
+(A + x n) • (B x n) ■/+(*)• fl(t))], 

(8.19) 



Pb {6,^:V,t) = 

(1 - F S ) ■ P B (9, i>, ip, t) + F S Q B (9, ^ <p, t) 

+2^Re [X M ((A_ x h) ■ (B x n) ■ |/_(t)| 2 
+(A + x h) ■ (B x h) ■ f + (t) ■ f*_(t))] . 

(8.20) 



9. Symmetries 

In this section we examine the symmetries of our differential rate formulae, starting from 
the simplest case, K + K~ in a P-wave, Eq. (3.8), but considering also the case where both 
P and S waves are included, Eq. (8.4). In the case of pure P-wave, one can readily spot 
that the probability densities in Eq. (3.8) are invariant to the following transformations: 

• A simultaneous rotation of the vectors A(t) and h 

• An inversion of the vector A(t) 

• Complex-conjugation of the vector A(t) 

The symmetry to simultaneous rotation of the vectors A(t) and n corresponds to the 
well-known freedom to choose a convenient basis in which to work. An example of an 
alternative basis is the helicity basis, which derives from the transversity basis by a cyclic 
permutation of the coordinate axis: xt = zh, etc. One can take the angles in Eq. (3.6) to 
be the polar and azimuthal angles in the helicity basis, but then one must transform A(t) 
accordingly, i.e, by permuting the elements of A(t) in the defining equation, Eq. (3.7). 
Then, Eq. (3.8) remains valid in the helicity basis. This rotational invariance implies that 
the choice of basis is irrelevant to the final result since the likelihood is invariant to the 
choice (though we do not rule out the possibility that the quality of the efficiency expansion, 
Eq. (4.2), may depend on the choice of basis, as pointed out in [14]). 

A more interesting symmetry is the symmetry that results from transforming A(i) 
to its complex conjugate. If we take, by convention, to be real and let 5» = arg(oii), 
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and 5± = arg(aj_), then as we will demonstrate below, this conjugation transformation is 
equivalent to the simultaneous transformation: 



(3 s tt/2 _ p a 



Ar 



-Ar 



5± — > it — 5± 
6\\ -4 27r - <5|| . 



(9.1) 



That is to say that the simultaneous transformation of these four variables is a symmetry 
of the likelihood because it transforms A(t) into its complex conjugate. Since for pure P 
wave state the combined transformation is a well-known symmetry, this observation may 
appear as a curiosity; however when both P and S wave states are included, we shall see 
that complex conjugation teaches us how to properly extend the symmetry. First, we show 
how the combined transformation accomplishes the claimed complex conjugation. 

1. Note from Eq. (3.4) that the combined transformation transforms E±(t) — > ±E±(t). 

2. Note also that the combined transformation transforms e~ 2t ^ s — > —e +2l ^ s and e +2ll3e —> 



3. Therefore, in Eq. (3.3), the terms in square brackets are transformed into their com- 
plex conjugates. 

4. Note that both cos 2/3 s and tl — t~h change sign under the transformation, so also the 
piece of Eq. (3.3) in the denominator, under the square root sign, is invariant under 
the combined transformation; since that piece is real we can say that it is anyway 
equal to its complex conjugate. 

5. The real quantity ao does not change under the combined transformation, but since 
it is real, it is anyway equal to Oq. 

6. The combined transformation transforms ay — > ajj. 

7. The combined transformation transforms ia± — > —ia* ± . 

8. Then looking at Eq. (3.7), one sees finally that the net effect of the combined trans- 
formation has been the complex conjugation of the vector A(t). 

Returning now to the full likelihood including both P and S wave states, Eq. (8.4), we 
can see that, here again, complex conjugation of the term 



leaves the probability density invariant (in a parameter space now enlarged to include 
Htf, and T^); however now, complex conjugation of the term g(/i), Eq. (8.5), implies that 
the transformation — > —T^ should also be carried out, in addition to the transforma- 
tion of (3 S , AT, Su, and 5±. Since negative values of are physically meaningless, this 
transformation is not an admissible symmetry. 



e -2i/3 s 




(9.2) 
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However we can find a symmetry transformation that carries one set of physically 
meaningful parameters into another. Such a symmetry is the transformation of the terms 
in Eq. (9.2) into their negative complex conjugate. This transformation is equivalent to the 
combined transformation already described, in addition to: 

o~ s — > 7T - S s 

(p-h) ( 9 - 3 ) 

The latter transformation carries us from a point on one side of the (ft mass peak to 
another point located symmetrically on the other side. This symmetry is useful when 
considering likelihood functions in which the dependence on \i is integrated out. If we 
integrate symmetrically about the <f> mass peak, we can consider the contribution to the 
integral coming from a slice in <f> mass on one hand and the a symmetrically-located slice 
in <p mass on the other hand. While the contribution of either slice is not invariant to 
the transformation above, the contribution of both slices certainly is, and the combined 
transformation: 

Ps -> vr/2 - p a 

Ar -> -Ar 

5± — > tt — 5± 

5 s ^ir-5 s (9.4) 

is again a symmetry of the integrated likelihood. We note, however, that this symmetry 
requires the symmetry of the nonrelativistic ^-propagator, Eq. (8.5), and applies to the 
likelihood integrated over a finite symmetric interval of integration. 

Symmetries of the likelihood function for B® — > J/ipfp? m t ne presence of S'-wave 
contribution for a fixed value of [i = m(K + K~) were discussed in a recent publication [15]. 
These formula can also used to fit for data falling within a narrow window in (M. Under 
those assumptions we can drop the <f> propagator from the expression in Eq. (9.2), absorb 
the Breit-Wigner terms into the amplitudes A(t), and write our model for the rates as 

y^Y s A(t) + e i5 > ^B(t) . (9.5) 

Then one can see that the transformation in which 5 S — > —5 S replaces 5 S — > it — 8 S in 
Eq. 9.4 accomplishes a complex conjugation of the terms in Eq. 9.5 and is a symmetry of 
the likelihood at fixed fi. 

In the more general case one can notice from Eqs. 8.19 and 8.20 that the probability 
densities integrated over [i are invariant to complex conjugation of both A(i) and the 
overlap integral X mu of Eq. 8.10. This can be accomplished with a more complicated 
adjustment of 5 S . With a nonrelativistic Breit Wigner the required transformation is 

6 S — > 25bw — 
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where 5 B w = arg (log ((^u - (1$ + ^<t>/ 2 )/ (Wo - M</> + ir V2))). The phase <5 B w reduces 
to <5bw = in the limit of an infinitesimally thin interval in fj,, and to 8bw = — vr/2 
in the limit of a finite symmetric interval. This demonstrates real differences in the two 
formulations, and underscores the need for caution when applying the formulae of Ref. [15] 
to a finite interval in fx = m{K + K~). 

10. Conclusion 

In this paper we have presented a compact formalism to easily access physical observables 
in the analysis of the decay B® — > J/ipcj). This formalism has practical applications, since 
complex vectors and their vector- algebraic operations can be easily implemented in high- 
level computer languages in order to model and generate such decays, but also because 
the symmetries of the formulae under operations such as rotation and complex conjugation 
are apparent and provide better physical insight into this complicated decay mode. The 
normalized probability densities can be used for the experimental extraction of physical 
parameters, in scenarios with no CP violation, with mixing-induced CP violation, or even 
with direct CP violation. In case of mixing induced CP violation, the effect of the S-wave 
contribution has also been included in the decay rate formulae. 
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